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1 . /ORE WORD  [ABSTRAa) 

A model  for  anomalous  self-diffusion  in  the  bcc-phase 
Group  IV  3 transition  metals,  Ti,  Zr  and  Mf  is  proposed.  The 
model  is  based  upon  the  identity  of  the  activation  complex  for 
vacancy  diffusion  in  the  bcc  lattice  and  the  smallest  unit 
(omega  embryo)  of  a hexagonal  phase,  the  omega  phase,  found  in 
these  metals  and  their  alloys  in  a metastable  state.  The  coef- 
ficient for  self-diffusion  is  calculated  as  a function  of  tem- 
perature by  means  of  a phenomenological  expression  for  the  free 
j energy  of  formation  of  an  omega  embryo,  and  by  means  of  a more 

fundamental  approach  based  on  Kikuchi's  Cluster  Variation  method. 
The  two  methods  of  calculation  for  the  self-diffusion  coefficient 


of  Zr  give  excellent  agreement  with  experiment,  the  observed 
curvature  of  the  Arrhenius  plot  being  reproduced  perfectly.  Further- 
more, diffusional  anomalies  observed  in  bcc-stabi 1 i zed  alloys  of  these 
metals  and,  in  particular,  the  increase  of  the  apparent  activation 
energy  with  solute  concentration  are  readily  explained  by  the  model. 

The  phenomenological  equations  derived  are  simple  to  use,  have  straight- 
forward  physical  content  and  may  be  employed  to  estimate  diffusion  coef- 


ficients in  a omalous  Group  IV  B transition  metals  and  alloys  when  little 
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3.  LIST  or  ILLUSTRAl IONS 


Figure  1:  The  con fi gurati ons  of  the  elementary  w embryo  and  the 

activated  complex  for  diffusion.  The  atom  at  0 exchanges 
place  with  a vacancy  at  V (a)  by  forming,  in  the  process, 
an  elementary  in  embryo  (b). 

Figure  2:  Comparison  of  the  diffusion  coefficient  predicted  by  the 

phenomenological  model  to  the  experimental  data  by  Federer 
and  Lundy  . 

Figure  3:  The  concentration  of  m embryos  (activated  complexes  for 

diffusion)  versus  the  absolute  temperature.  The  concen- 
tration of  ui  embryos  increases  with  a decrease  in  tempera- 

(15) 

ture  as  postulated  in  the  phenomenological  approach'  . 

Figure  4;  The  diffusion  coefficient  predicted  by  the  phenomenological 

(23) 

approach  is  compared  to  the  data  of  Murdock  and  Lundy  ' 

fpc 

for  self-diffusion  (a),  and  to  the  data  of  Santos  and  Dyment' 
for  diffusion  of  '*‘*Ti  in  Ti-Mn  alloys  of  composition:  (b)  9.7tMi 
(c)  13.3%Mn,  (d)  17.9.Mn  and  (e)  20.6';;.Mn. 
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4.  BODY  OF  REPORT 


4.  1.  Introduction 

The  objective  of  this  investigation  was  the  search  for  a 

theoretical  model  for  anomalous  self-diffusion  in  the  Group  IV  B 

transition  metals:  Titanium,  Zirconium  and  Hafnium.  These  three 

metals,  along  with  .-Uranium  and  the  rare  earths  6-Cerium  and 

n 2 31 

.-Praseodymium  exhibit  anomalous  diffusional  properties  ’ ’ . 

Normal  behavior  is  characterized  by  a diffusion  coefficient  D 
given  by 


D = Dp  exp(-0/P,T) 

2 

with  the  pre-exponential  factor  Dp  of  order  0.05  to  5.0  cm  /sec 

and  activation  energy  given  approximately  by  Q = 34  T^^  where  T^ 

is  the  absolute  temperature  of  melting  and  where  R and  T are  the 

gas  constant  and  absolute  temperature  respectively.  Dp  and  Q are 

approximately  temperature  independent  since  a plot  of  log  D vs. 

1/T  yields  a straight  line  over  several  decades  of  D.  By  contrast, 

the  anomalous  metals  are  characterized  by  significantly  curved 

Arrhenius  plots  which  have  been  interpreted  as  due  either  to  a 

temperature  variation  of  the  activation  energy  or  to  the  parti- 

(1  4) 

cipation  of  more  than  one  mechanism  in  bulk  diffusion  ’ . 

Numerous  attempts  have  been  made  to  explain  anomalous 

(1  4 5) 

diffusion  by  means  of  two  or  more  contributing  mechanism  ’ ’ . 

This  approdcfi  invariably  encounters  two  serious  difficulties: 


the  fact  that  anomalous  diffusion  is  confined  to  a very  specific 
qroup  of  metals  is  not  explained  nor  can  the  values  of  the  para- 
meters for  the  low  temperature  mechanisms  be  rationalized:  i.e., 
a pre-exponential  factor  which  is  several  orders  of  magnitude 

lower  than  that  of  normal  metals  and  an  activation  energy  at 

(1  4l 

least  half  as  large  as  the  high  temperature  one  ’ . 

The  approach  wh’rh  allows  for  a temperature  variation  of 
the  activation  en  -led  to  a few  investigators  to  be  a more 

(4  , 

sensible  one  :>uch  a temperature  dependence  of  Q is 

made  plausible  by  the  existence  of  a low  temperature  phase  transi- 
tion. Indeed,  all  anomalous  metals  present  allotropic  phase  tran- 
sitions; in  particular  Ti,  Zr,  Hf  and  Pr  transform  from  bcc  phase 
to  the  hexagonal  close  packed  structure  stable  at  low  temperature, 
and  the  stable  form  of  Uranium  at  room  temperature  is  orthorhombic, 
a distorted  hexagonal  close-packed  structure.  In  addition,  Ti , Zr 

and  Hf  transform  to  a metastable  hexagonal  (but  not  close-packed) 

(9l 

phase,  the  so-called  omega  phase,  upon  alloying'  ' or  under  pres- 
sure^^^^.  Furthermore,  Ti,  Zr,  Ce  and  Hf  are  Group  IV  B elements, 
the  transition  metal  counterparts  of  Group  IV  A elements,  the  light 
elements  of  which  are  convalently  bonded.  The  trigonal ly  bonded 
basal  plane  of  the  hexagonal  structure  has,  in  fact,  the  structure 
of  the  basal  plane  of  graphite  indicative  of  strong  covalent  compo- 
nent to  the  total  binding  energy. 
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Thus,  th('  bcc  structure  appears  to  be  a necessary  but  not 
sufficient  condition  for  the  diffusional  anomalies,  the  additional 
condition  required  being  that  the  bcc  structure  transform  to  a hexa- 
gonal phase  at  low  temperature.  Rather  than  seeking  a connection 
between  the  anomalies  and  the  hep  phase  itself^^\  we  attempted  to 
establish  the  connection  with  the  o phase. 

It  can  now  be  stated  unequivocally  that  the  objectives  of 
the  proposal  research  have  been  fully  met: 

1.  a model  which  explains  the  experimental  data  has  been  found 
and  gives  excellent  quantitative  fit, 

2.  subsequent  application  of  the  model  to  self-diffusion  in  Ti 
alloys^^^^  has  been  equally  successful, 

3.  our  model  is  claimed  by  other  investigators  to  be  the  most 

promising  one  to  explain  anomalous  diffusion  in  certain  bcc- 
(12) 

phase  rare  earths'  ' , 

4.  ours  appears  to  be  the  only  really  new  diffusion  model  proposed 

in  the  last  20  years.  As  such,  it  has  prompted  experimenta- 
lists^^ to  perform  new  measurements  on  the  anomalous  dif- 

fusers . 

A brief  summary  of  results  obtained  follows.  Details  will 

be  found  in  the  Ph.D  Dissertation  of  Dr.  J.  M.  Sanchez^^^^  and 

(15-17) 


in  relevant  publications 


4.?.  Diffraction  Model 


Our  first  task,  before  tackling  the  difficult  one  of  the 
diffusion  mechanism,  was  to  understand  better  the  oj  phase  trans- 
formation. We  therefore  attempted  to  explain  electron  and  neu- 
tron diffraction  patterns  from  bcc  Ti  and  Zr  alloys  transforming 
to  the  ^ phase.  It  was  found  essential  to  include  anharmonic 
effects  in  both  the  free  energy  and  the  diffraction  equations. 

In  this  way,  we  were  able  to  establish  a direct  connection  be- 
tween a Landau  expression  for  the  free  energy  of  the  transforming 
bcc  phase  and  the  diffraction  pattern.  Direct  Fourier  and  optical 

transforms  yielded  a plausible  picture  of  patterns  of  atomic  dis- 

(18) 

placements,  and  comparison  with  neutron  diffraction  data'  yiel- 
ded a plausible  temperature  dependence  of  the  third  order  anharmo- 
nic free  energy  term  B(T)  to  be  used  in  the  diffusion  model.  De- 


tails are  to  be  found  in  Ref.  16. 


8 


I 


4.3.  P h e n ome  no Jo^i  c d 1 Mo^o  1 

F roir  the  start,  we  assumed  a vacancy  mechanism,  i.e.,  a dif- 
fusing atom  exchanging  places  with  a vacant  substitutional  site,  'n 
so  doing,  the  atom  passes  through  an  unstable  saddle  point  configura- 
tion, momentarily  forming  an  "activated  complex"  which  spontaneously 
decays . 

The  key  to  the  explanation  of  anomalous  diffusion  consists  of 
recognizing  that,  in  bcc  structures,  the  required  activated  complex 
is  none  other  than  an  embryo,  that  is,  the  smallest  unit  or  clus- 
ter of  neighboring  atoms  having  the  essential  co-structure  characteris- 
tics. This  can  be  seen  in  Figure  1:  The  top  portion  (a)  shows  a 

(111)  plane  ABC  and  an  atom  at  0 occupying  the  normal  body-centered 
site  just  below  the  ABC  triangle.  For  this  atom  to  exchange  places 
with  a nearest-neighbor  vacancy  at  V,  it  must  first  pass  through  the 
"activated"  position  within  and  at  the  center  of  the  ABC  triangle. 

This  activated  configuration  is  shown  in  Figure  1 (b),  which  is  an 
. embryo  exhibiting  elementary  "collapse"  of  two  neighboring  (111) 
planes,  the  one  at  ABC  -and  the  one  at  0.  One  recognizes  in  Figure  1 
(b)  the  characteristic  trigonally  bonded  oj  structure'  . Once  the 
configuration  depicted  in  Figure  1 (b)  is  reached,  the  central  atom 
will  jump  to  a nearest  neighbor  vacancy,  if  available,  without  any 
extra  expenditure  in  energy  as  in  the  case  of  nonnal  diffusion.  The 
essential  difference  with  normal  diffusion  comes  about  when  the  re- 


I 

i 
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versiblo  work  of  tormdtion  of  the  activdted  state  is  calculated.  In 
the  case  of  systems  in  which  the  . structure  is  a low  temperature  me- 
tastable phase,  close  to  the  transition  temperature  the  free  energy 
of  formation  of  an  activated  complex  ' or  m embryo)  will  be  much 
smaller  than  in  the  case  of  normal  bcc  rrietals,  resulting  in  a higher 
concentration  of  activated  states  and  consequently  enhancing  diffusion 
at  low  temperatures.  On  the  other  hand,  in  the  limit  of  high  tempera- 
tures, random  thermal  vibrations  will  tend  to  destroy  the  u embryos 
and  the  diffusion  process  should  be  characterized  by  normal  values 
and  behavior  of  the  diffusion  coefficient. 

By  making  the  assumption  that  the  free  energy  change  in  forming 
unit  amount  of  is  described  by  a Landau  expansion  with  third-order 
coefficient  given  by 

B(T)  Const. /T  (1 ) 

it  is  possible  to  derive  the  following  expression  for  the  coefficient 
of  sel f-di ffusion 

D = Dq  exp(-Q/RT)  exp(AHJ^/Rn  (2) 

in  which  is  the  usual  preexponential  factor,  Q is  the  usual  activa- 
tion energy,  estimated  here  by  the  "normal"  rule  Q = 34  where  is  the 
melting  temperature.  is  the  motion  enthalpy  for  vacancy  diffusion 

and  is  the  bcc  ■)  transition  temperature. 
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In  the  particular  case  of  self-diffusion  in  bcc  Zr,  Q is 
70  kcal/mole,  and  is  estimated  to  be  about  890°K  by  extrapo- 
lation of  the  oj  transition  temperature  curve  vs.  composition  in 
(19) 

a Zr-Nb  alloy  . The  remaining  two  parameters  D„  and  .’.H  can 

0 m 

be  determined  by  fitting  Equation  (2)  to  the  data  of  Federer  and 
Lundy  yielding  = 0.21  cm^/sec  and  = 32.5  kcal/mole, 

both  parameters  falling  well  within  the  expected  range  of  values 
As  can  be  seen  from  Figure  2,  the  fit  is  as  perfect  as  the  data 
warrants.  Additional  details  are  given  in  Reference  15. 
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4.4.  C 1 u Vd  r i a t J on  Mode  1 

The  temperature  dependence  given  in  Lquation  (1)  was  merely 
postulated;  it  was  thus  desired  to  obtain  a more  fundamental  justifi- 
cation thereof.  To  that  end,  the  free  energy  of  the  bcc  phase  was 
written  as  a function  of  cluster  frequencies,  the  largest  cluster 
consisting  of  the  8-point  primitive  bcc  unit  cell,  each  point  being 
occupied  by  a 0,  •*!  or  -1  displacement  (on  a suitable  scale)  along  a 
given  [111]  direction.  It  was  thus  possible  to  simulate  activated 
complexes  or  other  inteniiedi ate  structure  between  perfect  bcc  and 
perfect  jj.  The  internal  energy  was  taken  as  a linear  combination  of 

cluster  frequencies  and  the  entropy  was  derived  by  the  Cluster 

(21 ) 

Variation  method  of  Kikuchi 

At  any  temperature,  equilibrium  states  are  obtained  by  mini- 
mizing the  free  energy  with  respect  to  cluster  frequencies.  This  is 
a formidable  task  since  there  are  3“'  = 6561  cluster  variables  (not  all 
distinct),  and  the  same  number  of  simultaneous  non-linear  algebraic 
equations  to  be  solved.  Dr.  Sanchez  was  able  to  do  this  by  modifying  the 
Natural  Iteration  algorithm  recently  developed  by  Kikuchi^^^^.  It  was 
then  possible  to  calculate  directly  the  temperature  dependence  of  acti- 
vated complexes,  as  shown  in  Figure  3 From  there,  it  was  possible,  as 
was  done  above,  to  obtain  the  temperature  dependence  of  the  diffusion 
coefficient  for  the  case  of  pure  Zr.  The  fit  to  the  data  by  this  method 
was  every  bit  as  good  as  that  shown  in  Fig.  2,  but  this  time  with  a solid 
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justification  for  all  temperature  dependences.  In  addition  to 
providing  new  insight  into  the  oi  phase  transformation , this  study 
provided  additional  returns;  (a)  it  was  apparently  the  first  ap- 
plication of  the  cluster  variation  method  (CVM)  to  displacive 
transitions,  perhaps  paving  the  way  for  future  use  of  the  CVM  to 
this  important  field  of  phase  transformations ; (b)  it  was  the  most 
ambitious  CVM  calculation  yet  undertaken,  in  terms  of  number  of 
variables  and  equations,  thereby  demonstrating  the  feasibility  of 
solving  very  complex  problems  by  the  algorithms  developed.  From 
the  strictly  computational  viewpoint,  the  experience  gained  has 
been  invaluable  for  our  new  ARO(D)-sponsored  project  which  has  just 
gotten  under  way. 

A complete  treatment  of  the  problem  will  be  found  in  the  Ph.D 
Dissertation  of  Dr.  Sanchez^^^^  and  in  forthcoming  publ ications . 
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.5.  Solvent  Diffosjon 

Group  IV  B trdnbition  metals  with  i- -s tabi  1 i z i nq  alloying  additions 
do  not  conform  to  well-established  rules  for  solvent  diffusion:  normally, 
if  an  alloying  addition  lowers  the  solidus  temperature  a corresponding 
lowering  of  the  activation  energy  for  diffusion  takes  place  as  expected 
from  the  relation  Q = 34  T^.  In  the  anomalous  metals,  however,  the 
effect  of  t- -stabi  1 i zing  elements  appears  to  be  the  exact  opposite:  al- 
loying additions  which  lower  the  solidus  temperature  result  in  an  increase 
in  the  activation  energy  for  diffusion.  For  instance,  solvent  diffusion 
in  Ti-rich  Ti-V  alloys  is  "anomalous'^  wi th  activation  energy  increasing 
with  V content,  while  V-rich  alloys  behave  normally.  Recent  experiments 
on  diffusion  of  ‘"‘Ti  in  Ti-Mn  and  Ti-Co  alloys^^^^  indicate  that  the 
curvature  of  the  Arrhenius  plot  decreases  as  the  alloying  element  content 
is  increased.  Also  if  the  data  is  analyzed  in  terms  of  a best  straight- 
line  fit,  both  effective  and  Q values  are  found  to  be  anomalously  low 
in  Mn  or  Co-lean  alloys,  to  become  progressively  more  normal  as  the  alloying 
element  content  is  increased. 

Such  behaviour,  which  defies  interpretation  in  terms  of  a normal 
diffusion  process,  is  readily  explained  in  terms  of  the  m-ernbryo  mechanism. 
Since  Mn(or  Co)is  a P-stabi  1 izi ng  element,  the  P >•  u transition  temperature 
decreases  rapidly  with  alloying  additions,  more  rapidly  in  fact  than 
the  solidus  temperature  This  means  that,  at  a given  temperature, 
increased  alloying  effectively  take  the  system  further  away  from  the 
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transition  so  that,  by  Figure  3,  the  concentration  of  di-einbryos  is 
expected  to  decrease,  resulting  in  increasingly  "noniial"  behaviour. 

The  phenomenological  model  presented  in  Section  3 is  particularly 
convenient  for  analyzing  these  effects  quantitatively.  By  consideration 
of  equation  (2),  we  may  define  an  effective  temperature-dependent  acti- 
vation energy 

T (c) 

Qgf^(T,c)  = 34  T,,,(c)[l  - . -f--  ] (3) 

where  . is  a constant  of  the  order  of  0.5,  and  where  both  solidus  and 
■ • <jj  transition  temperatures  and  are  concentration  dependent. 
Application  of  Equation  (3)  to  the  Ti-Mn  system  result  in  the  calculated 
curves  shown  in  Figures  4 a-e.  Not  only  is  the  fit  to  the  data  as  good 
as  can  be  expected,  but  the  resulting  values  of  the  fitted  parameters 
■dll  well  within  the  allowed  ranges,  while  the  curvature  of  the  Arrhenius 
plots  decreases  as  the  transition  temperature  decreases  with  concen- 
tration c.  Hence, It  can  be  concluded  that  all  anomalies  have  been  ex- 
plained Sdti s factori ly . 

Additional  details  are  given  in  the  Ph.D  Dissertation  of  Dr. 
Sanchez^'^^  and  in  forthcoming  publ ications ^ . 


4.6.  Conclusion 


Ihc  proposed  model  is  a simple  one.  witti  clear  (qualitative 
physical  ’ nterpretation  and  excellent  quantitative  fit  to  experi- 
mental data.  By  relatintj  the  diffusional  anomalies  of  Ti,  Zr  and 
Hf  metals  and  alloys  to  tlie  important  u displacive  transition,  we 
have  discovered  a new  diffusional  model,  surely  a significant  ac- 
complishment in  its  own  right.  Already,  our  model  is  cited  in 

(25  26 ) 

very  recent  review  articles'  ’ , and  is  described  as  "very  con- 

vincing" . 

From  a practical  point  of  view,  the  model  is  found  to  have 
predictive  value:  Fguations  (2)  or  (3)  are  sufficiently  simple,  and 
most  parameters  can  be  apprraximated  well  enough  to  enable  an  experi- 
menter to  predict  the  diffusional  properties  of  an  anomalous  system 
for  various  temperatures  and  concentrations,  without  the  necessity 
of  performing  costly  and  time-consuming  experiments.  For  example, 
resistence  to  high-temperature  steady-state  creep  in  Ti  alloys  can 
be  significantly  improved  by  adding  B-stabi 1 izing  elements,  and  by 
avoiding  hep-stabilizing  elements  such  as  A1 . The  analysis  can  be 
made  quantitative  through  the  use  of  Equation  (3). 

Much  wori  of  a fundamental  nature  remains  to  be  done.  Here 
are  a few  suggestions: 

a.  better  analytical  treatments  of  the  phase  transformation  are 
required,  either  in  the  framework  of  the  displacive  wave  method 
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with  annamonic  tonus  in  the  free  energy,  or  by  the  Cluster 
Variation  method  with  larger  clusters  than  the  ones  used  here. 

b.  the  best  quantitative  experimental  knowledge  of  the  , trans- 

fonnation  is  obtained  by  diffraction  methods;  X-ray,  electron, 

and  particularly  neutron  diffraction.  Much  more  work  is  needed 

to  collect  data  over  a range  of  concentrations  and  temperatures 

in  a variety  of  li,  Zr  and  Hf  alloys,  and  in  various  bcc-based 

(27 ) 

ordered  alloys  such  as  -CuZn' 

c.  Since  the  u phase  is  found  to  be  stable  in  the  pure  elements 

at  nigh  pressure,  experiments  of  both  diffraction  and  diffusion 
types  should  be  conducted  under  pressure  in  Ti  and  Zr  alloys 
as  a function  of  average  concentration  and  temperature. 

d.  other  anomalous  metals  such  as  ticc  Ce  should  be  investigated  in 
the  light  of  the  proposed  model  for  Ti,  Zr  and  Hf.  In  particular, 
it  would  be  interesting  to  find  out  whether  high-pressure  or  high- 
alloy  content  metastable  phases  with  structure  similar  to  <o  exist 
for  these  rare  earths.  Such  studios  have  not  been  conducted  pre- 
sumably because  these  elements  are  very  reactive  and  present  only 

a very  narrow  temperature  range  for  the  high-temperature  bcc  phases. 

e.  anomalous  diffusion  should  be  further  investigated,  in  both  element 
and  alloy,  by  such  techniques  as  the  isotope  effect,  as  is  currently 
being  performed  by  Professor  La/arus  and  co-workers.  Tarly  results 
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(14) 

appear  to  substatitiate  our  model  . 

In  the  past,  the  discovery  of  a new  theoretical  model  has 
often  suggested  renewed  experimental  activity  in  the  field,  while 
providing  useful  guid(!lines  for  practical  applications.  We  believe 
that  this  will  also  occur  in  the  present  instance  if  appropriate 
funding  is  available. 
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